Quantitative Test of SO (5) Symmetry in the Vortex State of Ndi^Ceo^CuO^ 
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By numerically solving models with competing superconducting and antiferromagnetic orders, we 
study the magnetic field dependence of the antiferromagnetic moment in both the weak and strong 
field regimes. Through a comparison with the neutron scattering results of Kang and Matsuura 
et al. on Ndi.gsCeo.isCuOi, we conclude that this system is close to a SO(5) symmetric critical 
point. We also make a quantitative prediction on increasing the upper critical field B C 2 and the 
superconducting transition temperature T c by applying an in-plane magnetic field. 
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PACS numbers: 74.20.-z, 74.20.De, 74.25.Op, 74.25. Ha 

A central question in the field of high T c superconduc- 
tivity concerns the nature of the competing state and its 
quantum phase transition to the superconducting (SC) 
state0, 0. The SO{5) theory^ considers the antifer- 
romagnetic (AF) state to be the dominant competing 
state, and describes its competition with the SC state 
in terms of an unified theory with approximate 5*0(5) 
symmetry. This theory predicts that the competing AF 
order reveals itself in the SC vortex core|3|,L|. This pre- 
diction has now been verified experimentally in a number 
of different compounds and by a variety of different ex- 
perimental techniques!! HBSIlEEIlElElEi- 
Experimental and theoretical progress on the AF vortex 
lattice is reviewed in recent articles jjj, |2j. 

While the experimental observation of enhanced AF 
correlation in the vortex state is encouraging, the limita- 
tion of low magnetic field could not convincingly estab- 
lish the AF state as the competing order. Only recently, 
it became possible to probe the magnetic structure of 
the AF vortex core with magnetic fields comparable to 
the upper critical field B c2 . Kang et al 5] and Matsuura 
et al\$ performed neutron scattering experiment on the 
electron-doped N d2- x Ce x CuO±-$ (NCCO) crystal un- 
der a magnetic field of 16T, beyond the upper critical 
field B C 2 = 6.2T of this material. The experiment finds 
a field induced AF scattering at the commensurate wave 
vector (7r,7r,0), where the AF moment scales approxi- 
mately linearly up to B C 2- The AF moment decreases 
with increasing field beyond B C 2- The experiment also 
finds field induced scattering at the (n, 0, 0) position. 
Scattering at this position has a different field depen- 
dence at high field, and may not be intrinsically related 
to the (7T, 7T, 0) scattering. We shall not discuss this scat- 
tering position further in this paper. 

The wealth of the experimental data on this system 
provides an opportunity to quantitatively test theoreti- 
cal models. In this work, we use the Landau-Ginzburg 
model of competing AF and SC order, including the Zee- 
man term. By quantitatively fitting the parameters to 
the experimentally measured field dependence of the AF 



moment, we could determine how close this system is to 
a 5*0(5) symmetric point. We show that the departure 
from the 50(5) symmetric parameters in general leads to 
departure from the linear dependence of the AF moment 
on the magnetic field, in the high field regime. There- 
fore, the experimentally observed linear dependence up 
to B C 2 directly implies that Ndi.^CeQ.i^CuOi is close to 
a 50(5) symmetric point. The magnetic field has a dual 
effect on the AF moment. On one hand, it creates AF 
moment in the vortex core, on the other hand, it reduces 
the AF moment by canting the spins uniformly towards 
to field direction. This explains the reduction of the AF 
moment when the magnetic field exceeds B C 2- We also 
show that the SC coherence length is self-consistently de- 
termined by the interaction of the SC and AF order pa- 
rameters, therefore, B C 2 is reduced in this system from 
its bare value. From this observation we predict that 
the value of B C 2 would be enhanced, if the AF order is 
suppressed by an in-plane magnetic field. The same rea- 
soning lead to an enhancement of T c in a system where 
AF and SC coexist uniformly. 

A minimal GL free energy that describes the competi- 
tion of AF and SC orders under magnetic field takes the 
form 



with 
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Here, we express the free energy and the order parame- 
ters in dimensionlcss form. is the standard GL free 
energy of the SC order parameter ip. T m is the free en- 
ergy of the AF order parameter m [15J . A is the vector 
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potential and b = (V x A) x a 2 is the dimensionless mag- 
netic field measured in units of cj>o/a 2 , where 4>o = hc/2e 
is the London flux quantum and a is the lattice constant. 

is the component of the AF moment perpendicular 
to the magnetic field, k is the GL parameter, fis is 
the Bohr magneton. The x term describes the Zeeman 
coupling of the AF moment with the external magnetic 
field. The interaction term u\i is the key term describ- 
ing the interaction between the AF and the SC order 
parameters. When u\ 2 > 1*21*1, there is a direct first 
order transition between the AF and the SC state. On 
the other hand, when u 2 2 < 1*2^1, the system undergoes 
two second order phase transitions, with an intermedi- 
ate phase where AF and SC coexist uniformly. In the 
special case of u 2 2 = u 2Ui, the potential can be scaled 
into a SO(5) symmetric form at the quantum critical 
point where r 2 j 'i/v® = fi/\/ui, and the AF and SC or- 
ders can be freely rotated into each other. Therefore, 
the relationship among the three quartic coefficients pro- 
vides a crucial test on the 5*0(5) symmetry of the model. 
On the other hand, by tuning the doping and the chem- 
ical potential, the quadratic coefficients can always be 
tuned to reach the quantum critical points. By quanti- 
tatively fitting the neutron scattering data in the vortex 
state of Nd\.^CeQ,\bCuO^, we can determine the value 
u i2/ ' \/u 2 u\, and consequently the nature of the quantum 
phase transition between AF and SC. 

We use the relaxation method^, 0] to find the nu- 
merical solution that minimizes the GL free energy 
We assume the system is an extreme type-II supercon- 
ductor, k > 1, so that the magnetic field is uniformly 
distributed over the sample. We choose the symmetrical 
gauge so that the vector potential A reads 



r x b 

A= - — 

The boundary conditions are then given by 

ijj(r + L x ±) = ip(r)e- i2nA * L * /a 

V>(r + L y y) = ^(r)e- l27r ^ i » /a 

m(r + L x x) = 77i(r) 

m{v + L y y) = m(r). 
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Here, L x and L y are the width and height of the magnetic 
unit cell that satisfy the condition (b z 4>o/a 2 )L x L y = (f)Q. 
A x and A y are the components of the vector potential 
A in the afe-plane. With the free energy Q and the 
boundary conditions (jlal) - (|4dll . we choose tp, tp* and m 
as independent variables. The relaxation iteration equa- 
tions read 



m (»+D = m W - S 2 ^ 
dm 



{n) 



(n) 
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FIG. 1: The plot of field dependence of AF moment for dif- 
feren parameters. The other parameters are p\ — p2 = a 2 , 
t"i = —1, f2 = —0.85, Ui = «2 = 1 and \ — 42.4. Here, the 
parameters are chosen such that the maximum SC order is 1 
and the SC coherence length at zero field equals the lattice 
constant a of our discrete model. 



where <5i and 82 are small positive numbers adjusted to 
optimize the convergence and n denotes the generation 
of the iteration. Starting from a proper initial state, this 
relaxation procedure quickly reaches convergence. We 
then repeat the relaxation procedure for different sizes of 
magnetic unit cell, or equivalently for different magnetic 
fields, to obtain the field dependence of the AF moment. 

When U12 = 0, the two order parameters decouple. 
The bare SC coherence length is given by £q = pi/\ri\, 
and the bare upper critical field is given by Bo = 
<^o/27t£q. However, when ui 2 > 0, the magnetic field 
induces AF moment in the vortex core regions, which 
renormalizes the coherence length and the upper criti- 
cal field to £ > £0 and B c2 < Bq. These values are 
self-consistently determined by solving the coupled GL 
equations in the vortex state. FigQ] shows the field de- 
pendence of the AF moment, for different values of U\2- 
We see that this value directly determines the curva- 
ture of the field dependence curve. For strong repulsion, 
u 2 2 > U2U1, the AF moment increases slowly at low field 
since the AF moment is strongly suppressed by the SC 
order in the bulk. At high field close to B c2 , the AF 
moment increases steeply since most of the SC order are 
destroyed. Therefore, the field dependence of the AF 
moment has a positive curvature in this regime. In con- 
trast, when u\ 2 < U2U1, the field dependence is reversed. 
The AF moment increases steeply for low field, and lev- 
els off at high field, leading to a negative curvature. In 
the SO(5) symmetric case where u\ 2 = u 2 u\, the AF 
moment scales linearly with B up to B C 2- 

The linear dependence was first predicted by Arovas 
et alQ in the low field regime, based on the argument 
that the AF moment scales with the number of vortices, 
which is linearly proportional to the field. Here we have 
shown that this linear dependence generally does not hold 
for fields close to B C 2, unless the system is SO(5) sym- 
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FIG. 2: The field dependence of the AF moment. The 
solid squares are obtained by subtracting the intensity of 
(0,0, 4.47t) measured at Hahn-Meitner-Institute (HMI)|2(| 
from the intensity of (tt, tt, 0) measured at National Insti- 
tute of Standards and Technology (j . The solid triangles 
are obtained subtracting the intensity of (0, 0, AAtt) mea- 
sured at HMl[20l from the intensity of (ir,ir,0) measured at 
HMI 6] . The parameters used in obtaining theoretical curve 
are pi — p2 = a 2 , n = — 1, ri = —0.85, wi = «2 = 1, 
U12 = 0.95 and x = 42.4. 

metric. Demler et al. |2| showed that the circulating 
currents around the vortices lead to a logarithmic cor- 
rection of the form B \n(B C 2/ B). Their result is derived 
in the weak field regime B <C B C 2. This behavior leads to 
an infinite slope at B = 0. The experimentally measured 
field dependence shown in Fig. 3e of reference []| does 
not have this feature. Zhang et al. studied the similar 
problem, in the B < B C 2 regime, mostly for the param- 
eters where u 2 2 < U1U2, and found the field dependence 
of the induced magnetic moment to have a negative cur- 
vature. This result is consistent with ours in the same 
parameter regime. 

For field higher than B C 2, the SC order is completely 
suppressed. Among those terms in Eq.Q, only the AF 
term T m and the magnetic field energy term survive. The 
problem becomes uniform and can be solved analytically 
to yield 



m(B) = y/\r 2 \/u 2 y/l- B*/B* 



where 



B„ 



(6) 



(7) 



is the critical field where the AF moment disappears. 

In FigEl we plot the fit of the experimental results 
of Matsuuraei <zZ@, EJ]. Both the intrinsic AF moment 
of NCCO and extrinsic FM moment of NCL2O3 impu- 
rity phase contribute to the (it, tt, 0) peak measured in 
HK0 zone [corresponding to field in (0, 0, 1) direction] [2]] 
. However, the extrinsic contribution of the NCI2O3 im- 
purity phase can be easily subtracted. Due to the cubic 
symmetry of NCI2O3 crystal, we can subtract the inte- 
grated intensity of (0, 0, 4.4-71") peak measured in HHL 



zone, which arises exclusively from Nd20%, from the nor- 
malized integrated intensity of (tt, it, 0) measured in HK0 
zone, to obtain the intrinsic NCCO contribution of the 
AF moment p3|. We then take the square root of the 
integrated intensity to obtain the moment. We obtain 
a reasonably well fit of the field dependence of AF mo- 
ment. For the field below B C 2, our theory predicts that 
the FM moment per vortex scales linearly with the field. 
It follows that the total FM moment scales quadratically 
since the number of vortices is proportional to the field. 
The experiment observed a linear relation. We believe 
the quadratic relation between FM moment and field is 
overwhelmed by the contribution of Nd moment, which 
is linearly polarized by the field. 

The quantitative fit to the field dependence of the AF 
moment determines U12/ ' y/uiu\ = 0.95. Therefore, the 
quartic terms are approximately 50(5) symmetric. The 
quadratic terms r% = — 1 and T2 = —0.85 deviate more 
from the SO(5) symmetric value. This is to be expected, 
since the system is close to optimal doping where the 
SC state is the stable ground state. Since n depends 
on the chemical potential, it can always be tune to r 2 at 
the quantum phase transition point where both quadratic 
and quartic terms are approximately 5*0(5) symmetric. 
The experimental result of the field dependence of AF 
moment beyond B C 2 can be fitted very well by the sim- 
ple formula Eq.© with a critical field B s ~ 20. 5T. En- 
ergetically, this critical field is determined by the com- 
petition between the Zeeman energy and the conden- 
sation energy of the AF moment. From the fact that 
the observed AF moment at B C 2 is about m — 0.05/is 
per sitefEj], we estimate that the Zeeman energy at B s is 
about B s m w 0.05met>. It follows that the AF coupling 
J af is about 20met> , if the AF condensation energy can 
be written as Eaf = Jaf^ti 2 - 

The detailed agreement between our model and the 
experiment obtained above allows us to make a striking 
quantitative prediction. As we mentioned previously, the 
upper critical field B C 2 is strongly reduced by the repul- 
sion between AF and SC. Using the parameters obtained 
above, we estimate the bare value of the SC coherence 
value to be t^Q = pi/|ri|, and the bare upper critical field 
is given by B = c/) \r 1 \/2irp 1 w 23T. If the AF order 
could be suppressed without disturbing the SC order, the 
SC order would be enhanced. Consequently, the upper 
critical field B C 2 would increase too. Our idea is to apply 
a Zeeman field in the afe-plane. The Zeeman field sup- 
press the AF moment through the canting term \- After 
including the Zeeman field, the total field is 



b = b 7 i + 6iix. 



(8) 



The dimcnsionless upper critical field b C 2 is then given by 
the solution of the coupled equations: 



b 2 c2 -H 



It I - u 12 m 2 
2irpi/a 2 
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FIG. 3: The angle dependence of the z-component of the 
upper critical field, given by the solution to coupled equations 
l|9afl and l|9b|l . The parameters are the same as those used in 
fitting of the experimental data in Fig|5| 
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The first equation gives the z-component of 6 c2 as a func- 
tion of to while the second one gives the AF moment m at 
b C 2- For small small fry, the increase of the z-component 
of b C 2 is quadratic in 611 



b*Ah)- b c2 (Q) 



X ui2bu 



0(bf,) 



u 2 \ri\ - u 12 \r 2 \ - XUi 2 (b z c2 (0)y 



(10) 



where b z c2 (0) is the upper critical field if b\\ = 0, given by 
the solution of self-consistent equation 



&(0) = 



u 2 \ri\ - u 12 \r 2 \ + XMi2(bg 2 (0)) 
2irii2Pi/a 2 



(11) 



This effect can be used to quantitatively measure the 
competition between the AF and the SC order parame- 
ters. FigOU plots the z-component of the upper critical 
field for NCCO system studied in previous paragraphs, 
as a function of the canting angle 9 = tan _1 (&* 2 /fr||)- ft 
is interesting to note that a canting of 27° increases the 
B* 2 by about 30%. fn the above estimate we neglected 
the effect of the ab plane field on the SC. If this effect is 
taken into account, the actual increase of B z 2 would be 
smaller. 

Similarly, for a system with uniformly co-existing AF 
order and SC order at zero field, the suppression of the 
AF ordering due to the Zeeman field in the ab plane 
would increase the transition temperature T c itself, ac- 
cording to this formula: 



r c (6u)-r c (o) 

T c (0) 



XUubf, 



(12) 



Here, we assume r 2 (T) = r 2 (l - T/T%), n(T) = r x (l - 
^ with T^r > T°. We also assume that the Zeeman 



field 6|| is small enough so that the AF moment m is 
nonzero at T c . Again, the enhancement scales quadrati- 
cally with the Zeeman field. 

In summary, we have studied the competition of AF 
and SC orders in High-Tc superconductors by solving the 
GL free energy. We showed that the curvature of the AF 
moment versus B field plot directly determines the de- 
parture from the 5*0(5) condition in the quartic terms of 
the GL functional. Reasonable agreement with the recent 
neutron scattering experiment on Ndx^Ce^.x^CuOi 
shows that this system is close to a 5*0(5) symmet- 
ric quantum critical point. New experiments are pre- 
dicted to increase B c2 and T c by applying an in-plane 
magnetic field. The GL model with competing AF and 
SC order parameters can be derived from the projected 
SO (5) model on the lattice p^. which can in turn be 
derived from the microscopic t-J model by a contractor- 
renormalization-group algorithm |24j. The quantitative 
agreement between these approximately 50(5) symmet- 
ric models and the experiment shows a promising direc- 
tion towards a full microscopic theory of high Tc super- 
conductivity. 
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